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ABSTRACT 

We analyse five observations of the X-ray binary GX 339-4 by the soft 7-ray OSSE 
detector on board CGRO simultaneous with either Ginga or RXTE observations. 
The source was bright during four of them, with the luminosity of L ~ 10 37 erg s _1 
and the spectrum typical for hard states of accreting black holes, and it was in an 
off state during the fifth one, with L ~ 10 35 erg s _1 . Our broad-band spectral fits 
show the mean electron energy of electrons in the Comptonizing plasma decreasing 
with increasing luminosity within the hard (bright) state. For the observation with the 
best statistics at soft 7-rays, ~ 1/4 of energy in the Comptonizing plasma is probably 
carried by non-thermal electrons. Then, considering the efficiency of Comptonized 
hybrid synchrotron emission allows us to obtain an upper limit on the strength of 
the magnetic field in the X-ray source. Furthermore, this synchrotron emission is 
capable of producing the optical spectrum observed in an optically-high state of GX 
339-4. In the off state, the hard X-ray spectrum is consistent with being dominated 
by bremsstrahlung. The unusually strong Fe Ka line observed by the RXTE/PGA 
during that state is found not to be intrinsic to the source but to originate mostly in 
the Galactic diffuse emission. 

Key words: accretion, accretion discs - binaries: general - stars: individual: GX 
339-4 - gamma-rays: observations - gamma-rays: theory - X-rays: stars. 



1 INTRODUCTION 

The Galactic X-ray source GX 339-4 was discovered by the 
OS 0-7 satellite (Markert et al. 1973). It exhibits state tran- 
sitions typical for accreting black holes in binaries; it has 
been observed in the quiescent, off, hard (low), soft (high) 
and very high X-ray states. Its X-ray spectral analyses have 
been done, e.g., by Ueda, Ebisawa & Done (1994, hereafter 
U94), Zdziarski et al. (1998, hereafter Z98), Wilms et al. 
(1999) and Kong et al. (2000). In the soft 7-ray domain, 
OSSE observations have been studied by Grabelsky (1995), 
Z98, and Smith et al. (1999), and SIGMA results have been 
presented by Bouchet et al. (1993) and Trudolyubov et al. 
(1998). X-ray variability studies have been done, e.g., by 
Maejima et al. (1994), Smith & Liang (1999), Nowak, Wilms 
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& Dove (1999) and Revnivtsev, Gilfanov & Churazov (2001, 
hereafter R01). 

The source has been also observed in the radio (Fender 
et al. 1999 and references therein) and an extended emission 
from a jet-like elongated structure was reported (Fender et 
al. 1997). Optical observations (e.g., Grindlay 1979; Motch, 
Ilovaisky & Chevalier 1982; Motch et al. 1983, hereafter 
M83; Soria et al. 1999) have not been able to identify the 
donor star (see also a discussion in Z98). In fact, the ob- 
served optical radiation appears to be mostly produced in 
an accretion flow and/or an outflow. This is supported by 
the decline of the optical flux during a transition from a hard 
to a soft state (M83; Makishima et al. 1986). Interestingly, 
this type of transition is also accompanied by a decline of 
the radio flux (Corbel et al. 2000). 

The presence of an accreting black hole in the system is 
also suggested by its Fourier power spectra having the same 
characteristics at high frequencies as other black-hole candi- 
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dates, but not neutron-star systems (Sunyaev & Revnivtsev 
2000). Also, its broad-band energy spectra are remarkably 
similar to those of a well established black hole binary, Cyg 
X-l (e.g., Gierlihski et al. 1997, 1999). However, no reliable 
constraints on the mass of the compact object exist (see a 
discussion and references in Z98). 

X-ray and 7-spectra of black holes in their hard states 
are successfully described by models with thermal Comp- 
tonization of soft photons in optically thin plasmas of a 
temperature, T ~ 10 9 K, and reflection of the Comptoniza- 
tion spectrum from an optically thick, cold (T ~ 10 6 K) 
medium, presumably an accretion disc (Sunyaev & Triimper 
1979, Zdziarski et al. 1995, Gierlihski et al. 1997). The disc is 
also likely to serve as the source of soft photons undergoing 
Comptonization. This view is supported by the correlation 
between the spectral index of the Comptonization spectrum 
and the amount of reflection observed in GX 339-4 (U94, 
R01) and in other accreting black holes (Zdziarski, Lubihski 
& Smith 1999; Gilfanov, Churazov & Revnivtsev 2000). 

The role of magnetic helds remains an open question. 
They are probably involved in transfer and generation of 
energy either via viscosity mechanisms or reconnection pro- 
cesses. As a consequence, synchrotron radiation should be 
present and it may serve as a source of soft photons under- 
going further Comptonization. Though for thermal plasmas 
this process appears to be typically negligible (Wardzinski 
& Zdziarski 2000, hereafter WZ00), the presence of a weak 
non-thermal tail strongly amplifies synchrotron emission 
(Wardzinski & Zdziarski 2001, hereafter WZ01). Then, the 
synchrotron process can become a very effective source of 
soft photons, comparable to or stronger than the cold disc 
emission. 

Consequently, detection of a non-thermal electron tail 
via a high-energy tail in a Comptonization spectrum can put 
constraints on the strength of the magnetic field. In particu- 
lar, observations of Cygnus X-l in its hard state suggest the 
presence of such a tail (McConnell et al. 2000, 2002). For 
this case, WZ01 obtained an upper limit on the magnetic 
field strength in the accretion flow and found constraints on 
the geometry. 

Here, we study those (and some other) issues in the 
context of GX 339-4. We analyze five Ginga and RXTE ob- 
servations of GX 339-4 which were simultaneous (or nearly 
simultaneous) with 7-ray observations by the OSSE detec- 
tor onboard CGRO. This provides us with data in the useful 
range from a few keV to a few hundred keV. The data re- 
duction procedures and models we use are described in Sec- 
tions U and ^, respectively. We present the results of spectral 
analyses in Sections ^ and ^. Consequences of the possible 
presence of a non-thermal electrons for accretion flow mod- 
els are discussed in Section ^j. We investigate the connection 
between the X-ray and optical data in the framework of 
synchrotron self-Compton models in Section [?], while X-ray 
variability is examined in Section ^. Our results are summa- 
rized in Section ^. 

2 SELECTION AND REDUCTION OF 
OBSERVATIONAL DATA 

The data used in this work are listed in Table 1 . The first two 
data sets are identical with those of Z98. They consist of two 



Ginga data sets from 1991 September 11 and 12, with the 
usable energy range of 1.2-30 keV, and simultaneous OSSE 
data in the range from 50 keV up to several hundred keV. 
A 1 per cent systematic error is added in quadrature to the 
statistical error in each Ginga channel (as in U94). All the 
OSSE data used here include energy-dependent systematic 
errors estimated from the uncertainties in the low-energy 
calibration and response of the detectors using both in-orbit 
and prelaunch calibration data. They decrease from ~ 3 per 
cent at 50 keV to ~ 0.3 per cent at > 150 keV. 

The third data set consists of a relatively short RXTE 
observation on 1996 July 26, which was performed shortly 
after a two-week long OSSE observation of 1996 July 9-23. 
For these RXTE data (and other ones used here), we used 
data from the units and 1 of the PCA detector, and used 
only layer 1, which provides the highest signal-to-noise ratio. 
This choice is favoured by Wilson & Done (2001) based on 
their study of Crab spectra from the PCA. Also following 
them, we added a systematic error of 0.5 per cent to the 
PCA data. The RXTE data were reduced with the lhea- 
SOFT package v. 5.0.4, response matrix v. 7.10. Standard 
data selectrion criteria were applied, namely the Earth ele- 
vation angle > 10°, pointing offset < 0.02°, the time since 
the peak of the last SAA passage greater than 30 min and 
the electron contamination < 0.1. We used the Sky_VLE 
background model for bright sources dated 1999 September. 
We constrain the energy range to a conservative range of 
3-30 keV. We also used data from the clusters A and B of 
the HEXTE detector in the 30-100 keV range. 

The fourth data set consists of an RXTE observation on 
1997 February 14 and an OSSE observation on 1997 Febru- 
ary 11-19. This data set has the best statistics among all 
ones used here. 

Finally, the fifth data set includes an RXTE observation 
on 1999 July 7 and a two-week long OSSE observation of 
1999 July 6-13. In the PCA energy range, the source was 
by two orders of magnitude fainter than in the other data 
sets. Thus, in order to obtain the best signal to noise ratio, 
in the PCA data analysis we used the layer-1 data from all 
the detector units available during this observation, i.e., 0, 
2 and 3. Here we employed the background model for faint 
sources from 2002 February. The HEXTE data were not 
used because of very large statistical errors after background 
subtraction. 



3 SPECTRAL MODELS 

Fits below are performed with the xspec package (Arnaud 
1996) v. 11. The confidence ranges of each model param- 
eter are given for a 90 per cent confidence interval, i.e., 
Ax 2 = 2.71 (e.g. Press et al. 1992). On the other hand, 
the plotted vertical error bars are 1-a, the upper limits, 2-tr, 
and the plotted spectral data are rebinned for clarity of the 
display. Spectra from various detectors are allowed to have 
free relative normalization. 

In most cases, the intrinsic continuum is modelled by 
Comptonization of blackbody photons in a hybrid plasma 
with the COMPPS model of Poutanen & Svensson (1996). 
Compton reflection of the continuum from an ionized, opti- 
cally thick medium (presumably a cold accretion disc) is also 
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Table 1. A summary of the observational data. The background-subtracted count rates for the LAC, PCA, HEXTE and OSSE 
instruments arc for the 1.2-29, 3-30, 30-100 and 50-1000 kcV energy ranges, respectively. 



No. 


Satellite 


Date 


Observation ID/ 


Instrument 


Exposures 


Count rate 


References 








Viewing Period 






[s 




Is" 1 ] 




1 


Ginga 


1991 Sept. 11 


91091109421 


LAC 


2.0 


X 


10 3 


1720 


U94; Z98 




CGRO 


1991 Sept. 11 


VP 9 


OSSE 


4.8 


X 


10 3 


15 


Grabelsky et al. (1995) 


2 


Ginga 


1991 Sept. 12 


91091109421 


LAC 


3.8 


X 


10 2 


1700 


U94; Z98 






i noi g«t-,+ i o 
iyyi oept. iz 


\rT> n 
Vr y 




2.5 


X 


10 :1 


10 


oraDeisRy et al. (lyyoj 


3 


RXTE 


1996 Jul. 26 


P10420-01-01-00 


PCA 


5.2 


X 


10 3 


271 


Bottcher et al. (1998); 










HEXTE A 


1.7 


X 


10 3 


16 


Smith & Liang (1999); 










HEXTE B 


1.7 


X 


10 3 


13 






CGRO 


1996 Jul. 9-23 


VP 524 


OSSE 


5.4 


X 


10 r> 


6.5 


Smith et al. (1999) 


4 


RXTE 


1997 Feb. 14 


P20183-01-02-00, 01 


PCA 


2.1 


X 


10 4 


217 


R01 










HEXTE A 


6.5 


X 


10 3 


13 












HEXTE B 


3.6 


X 


10 3 


11 






CGRO 


1997 Feb. 11-19 


VP 614.5 


OSSE 


3.4 


X 


10 r> 


7.0 




5 


RXTE 


1999 Jul. 7 


P40108-03-01-00 


PCA 


1.2 


X 


10 4 


1.5 


Feng et al. (2001) 




CGRO 


1999 Jul. 6-13 


VP 821 


OSSE 


1.1 


X 


10 s 


0.2 





a For HEXTE, the detector live time (i.e. dead-time corrected) is given. 



taken into account using the inclination-dependent Green's 
functions of Magdziarz & Zdziarski (1995). 

In the case of a purely thermal plasma, the Comp- 
tonization spectrum is characterized by the Thomson op- 
tical depth, r, the electron temperature (which we express 
in energy units, kT, where k is Boltzmann constant), and 
the temperature of the blackbody seed photons, fcT see d. We 
also consider a hybrid electron distribution, which consists 
of a Maxwellian up to a electron Lorentz factor, 7 n th, and a 
power-law of an energy index, p, above it. Given these pa- 
rameters, we calculate the ratio, S, of the energy in the non- 
thermal electrons to that in the thermal ones. In order to 
quantify the hardness of a given Comptonization spectrum, 
we also calculate the 4-20 keV power-law photon index, V, 
which we obtain by fitting data simulated with the model of 
(only) Comptonization. 

The reflection is characterized by the solid angle sub- 
tended by the reflecting medium, Q, its ionization parame- 
ter, £, and the Fe abundance relative to the solar one (Anders 
& Ebihara 1982), Af b . Here, £ = Li on /nR 2 , L- 10 n is defined 
as the 5 eV-20 keV luminosity in a power law spectrum and 
n is the density of the reflector located at distance R from 
the illuminating source. We use here the ionization model 
by Done et al. (1992), which is not applicable to very highly 
ionized Fe (Ballantyne, Ross & Fabian 2001), but sufficiently 
accurate at low/moderate ionization, such as that found in 
our data. The reflector temperature is kept at 10 6 K, which 
is about the highest temperature consistent with the model 
of ionization equilibrium of Done et al. (1992). Compton 
reflection is accompanied by a fluorescent Fe Ka line at a 
rest-frame energy, Ef , and an equivalent width, Wfc- Both 
the line and reflection are relativistically broadened in an ac- 
cretion disc in the Schwarzschild metric with the emissivity 
equal to that of a thin disc, extending from an outer radius, 
r ut (assumed to equal 10 3 ), to an inner one, r ln > 6, where 
r is in the units of GM/c 2 (Fabian et al. 1989). We note 
here that the values of kT, j nt h, and p can be constrained 



mostly by the OSSE data whereas r (for given kT) and the 
reflection/line parameters are determined mostly by either 
Ginga or RXTE data. 

We assume the column density to the source of Nh = 
6 x 10 21 cm' 1 . This value was derived from E(B — V) mea- 
surements and is consistent with results of other methods, 
see Z98 and references therein. We assume the distance to 
the source of d = 4 kpc (which agrees with both kinematic 
and extinction measurements, see Z98). In calculations of 
reflection and relativistic broadening, we assume an incli- 
nation of i — 45°, and we also assume M = IOMq, but 
we stress that currently there are no constraints on these 
quantities. 



4 HARD STATE OBSERVATIONS 

The source flux (extrapolated to the range of 1-1000 keV) 
in the data sets from 1 to 4 is F ~ 10 -9 erg cm -2 s _1 , which 
corresponds to the luminosity of L ~ (1-3) X 10 37 erg s _1 
and ~ 1 per cent of the Eddington luminosity (at d = 4 kpc 
and M = IOMq assumed here). 

4.1 The Ginga-OSSE data 

Results of thermal-Comptonization fits to the Ginga-OSSE 
data sets 1 and 2 are given in rows 1 and 2 of Table ^, respec- 
tively. The former has much better statistics and we discuss 
it in more detail, but results for the set 2 are very similar. 
As it was shown by Z98 (see their section 5.1), spectral fits 
to the Gint/a-OSSE data do not constrain the geometry of 
the Comptonizing plasma and the seed-photon source. Thus, 
we assume here a simple geometry of a spherical source with 
seed photons emitted uniformly. We also account for the soft 
excess present in the data by a blackbody component at the 
same temperature as that of the seed photons. 

We obtain a very good fit to the data set 1, with 
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Figure 1. (a) Broad- band spectra of GX 339-4. The data sets 
1, 3 and 4 (with decreasing flux) are shown in the green, red and 
blue error bars, respectively. The black solid curves give the best 
thermal fits to the data, and the dashed curve gives the best non- 
thermal fit to the set 4. (b) The spectral components of the model 
of the set 1 shown separately. Thermal Comptonization, the soft 
excess, Compton reflection and the Fe K line are shown by the 
long-dashed blue, dot-dashed red, short-dashed green and dotted 
magenta curves, respectively. The solid black line gives the total 
spectrum. Absorption by the interstellar medium is seen below a 
few keV. 



X 2 ' jv — 46/81. The value of x 2 l v much smaller than unity 
results from relatively large systematic errors added to the 
Ginga data, see Section ti. The fit components, as well as 
comparison of this spectrum with spectra from data set 3 
and 4 (analyzed in Section [T^ below) are shown in Fig. 
|l[ Fig. g shows the quality of our fits to the data sets 1- 
4. The parameters of the fit to data set 1 are very similar 
to those obtained by Z98 and indicative of a typical hard- 
state spectrum of black-hole binaries. The main difference 
with respect to the results of Z98 in their table 2 is the 
plasma temperature, which we fit here at kT = 46^4 keV, 
whereas Z98 obtained kT = 57lg keV. This difference is 
due to the model of Comptonization used here being more 
accurate than that used by Z98. 

The fitted relative normalization of OSSE to Ginga is 
close to unity, O.QltJJ;^ and 0.90^0 \ for the sets 1 and 2, 
respectively. Hereafter, we give the luminosity normalized to 
the OSSE data. Such a choice is appropriate for investigating 
a c onnection between the luminosity and kT (see Section 
4.3) since both are determined mostly by the OSSE data. 

In non-homogenous geometries (such as a corona above 
a disc), an anisotropy break appears at the peak energy of 
the second-order scattering for a high enough kT, with defi- 
ciency of photons below the break with respect to the extrap- 



olation of the main power law (e.g., Poutanen & Svensson 
1996). In general, the form or the absence of this break may 
be then used to constrain the geometry. However, the energy 
range where it may appear in the case of GX 339-4 is similar 
to that of the soft excess. Consequently, fits assuming differ- 
ent geometries and including the soft excess give very similar 
statistics [see also discussions in Z98 and Nowak, Wilms & 
Dove (2002)]. The presence of the soft excess itself is very 
significant statistically, e.g., Ax 2 = +29 for a fit to the set 
1 in spherical geometry and without the soft excess. 

We also note that our adopted simple model of the soft 
excess as a blackbody is probably not realistic. In particu- 
lar, the fitted fcT aeo d ~ 0.3 keV is about twice higher than 
that of Cyg X-l in the hard state (Ebisawa et al. 1996). In 
Cyg X-l, the soft excess is consistent with the presence of 
both blackbody and its Comptonization in a plasma cloud 
separate from the main Comptoning source (Di Salvo et al. 
2001; Frontera et al. 2001a). This results in the onset of the 
excess at ~ 3 keV (similar to GX 339-4) , but with the excess 
~ 1-3 keV emission being mostly due to Comptonization, 
not blackbody. Consequently, we also have fitted the soft 
excess by a second thermal Comptonization component (in 
addition to the main Comptonization component responsi- 
ble for the emission at > 3 keV) . We find the excess can be 
fitted by emission of either an optically-thick plasma with 
r = 6.1 and kT = 1.4 keV (using COMPTT, Titarchuk 1994), 
or optically-thin one with r = 0.43, kT = 20 keV (using 
COMPPS). Both models yield x 2 l v = 45/79, and are anal- 
ogous to the models of the soft excess in Cyg X-l by Di 
Salvo et al. (2001) and Frontera et al. (2001a), respectively. 
We conclude that the data do not allow us to distinguish 
between different soft excess models. 

As far as the high-energy end of the measured spectrum 
is concerned, an important issue is the possible presence of a 
high-energy tail in the electron distribution. We have found 
that none of the OSSE data considered in this work have 
enough statistics at high energies to constrain the electron 
power law index, p. We have thus chosen p = 4, which is 
close to the values measured in Cyg X-l, see Gierlinski et 
al. (1999), and McConnell et al. (2000, 2002), and it was 
used in the theoretical study of WZ01. 

We find that the Ginga-OSSE data are consistent with 
the presence of a high-energy electron tail, but they do not 
constrain its amplitude. E.g., the non-thermal fraction, S, for 
the data set 1 is 0.031q'q3 ■ The presence of a tail is indicated 
by the result of Z98, who obtained A^ 2 = —4 when allowing 
for the presence of non-thermal electrons in their fit to the 
OSSE data from 1991 September 5-12. However, an alter- 
native explanation of the departure of the spectrum for the 
thermal one can be a distribution of the plasma parameters, 
in particular kT, either in time or in space. 



4.2 The RXTE-OSSE data 

Since the usable PCA data extend down only to 3 keV (com- 
pared to 1.2 keV for the Ginga data), the soft excess cannot 
be constrained by them at all. Thus, we neglect its pres- 
ence for all RXTE-OSSE data. Then, the thermal model 
applied to the data set 4 (fit 4a in Table ^) yields an ac- 
ceptable x 2 l v ~ 143/148. The normalization of the models 
to the HEXTE A and B clusters relative to the PCA data 
is 0.88tp;o3 and 0.85l°' £3, respectively. On the other hand, 
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Table 2. Main parameters of the fits in Section ^. The numbers correspond to the data sets in Table H and a and b denote the thermal 
and hybrid fits, respectively, to the set 4. kT, L(l-1000 keV), £, -Efcj an d Wp a are is units of keV, 10 37 erg s — x , erg s — 1 cm, keV and 
eV, respectively. 



Fit 



kT 



Q/2n 



E Fc 



1 
2 
3 

4a 
4b 



46 
43 
58 
76 



+6 
4 
+ 8 
(i 
+3 
3 
+ . r . 



16 



+ 11 



2.5 
2.3 
0.88 
1.2 
1.2 



2.5 
2.7 
1.9 
1.5 
2.2 



0.3 
0.3 
+0.4 
0.3 
+0.1 
0.1 
+0.1 
0.1 
+0.4 



1.75 
1.73 
1.79 
1.79 
1.76 



0.36 
0.21 
0.40 
0.39 
0.27 



+0.07 
0.07 
+0.09 
0.09 
+0.06 
0.06 
+0.04 
0.03 
+0.05 
0.06 



130 

340 
110 
90 

330 



+ 110 
100 
+ 1050 
-260 
+70 
40 
+30 
30 
+ 460 
170 



2.9 
1 

3.6 
2.8 
2.3 



1.7 
0.9 
+3.6 
-0.9 
1.1 
1.8 
(1.8 
0.6 
0.8 
0.7 



0.25 
0.27 
0.21 
0.16 
0.33 



0.02 
0.03 
+0.03 
0.03 
+0.06 
0.10 
+0.05 
0.15 
+0.04 
0.06 



6.58 
6.58 
6.28 
6.31 
6.23 



0.32 
0.32 
+0.30 
0.32 
+0.14 
0.14 
+0.12 
0.13 
+0.11 
0.10 



70 
61 

88 
59 
105 



4M 

34 
+34 
33 
+ 17 
-18 
11 
12 
+ 11 
13 



36] 
185" 



30 



179 
6 +65 

°° — 27 

6 +io 



46/81 
63/81 
159/130 
143/148 
133/147 



X 



5 


-5 
5 



-5 
5 





-5 
5 
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Figure 2. (a) The contribution to the total x 2 from data bins multiplied by the sign of (data- model) and (b) data-to-model ratios for 
fits presented in Table ^. Red lines represent Ginga or PCA data, while green, magenta and blue lines correspond to the HEXTE A, 
HEXTE B and OSSE data, respectively. 



the relative normalization of models to the OSSE and PCA 
data was almost unity, 0.98]too3> m spite of the observing 
time being larger for the OSSE. The results of our fits in 
the X-ray regime are similar to those of Wilms et al. (1999), 
who analyzed an RXTE observation that took place three 
days after ours. 

Our data set 4 has the best quality of all those con- 
sidered here, and allows us to study the presence of a non- 
thermal electron tail in the source. We have found it im- 
proves the fit (4b in Table §) very significantly, to x 2 /v = 
133/147, with the probability of the fit improvement be- 



ing by chance of only 10 -3 (using the F-test, Bevington & 
Robinson 1992) . The parameters of the hybrid electron dis- 
tribution are 7 nt h = 1.59^014 and kT = 46lg X keV, which 
corresponds to 5 = 0.34^o'i7 (5 increases with decreasing 
either kT or 7, lt h). The normalizations of the models to the 
HEXTE A and B clusters and to the OSSE detector, relative 
to the PCA data, were 0.86±o.o3, 0.84t^o3 and 0.92±° ° t, re- 
spectively. 



We caution, however, that the OSSE observation was 
much longer than the RXTE one, and thus it is possible that 
the modification of the shape of the OSSE spectrum was due 
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not to non-thermal electrons, but resulted from variations of 
the plasma parameters in time or in space. We have checked 
these possibilities. First, we allowed for different values of 
r in thermal fits to the RXTE and OSSE data. This, how- 
ever, did not change the fit quality, y 2 /v = 142/147, and 
the difference between the two values of r was only 0.1. We 
then fit two COMPPS models allowing for different values of 
r and kT and keeping all other model parameters the same. 
This significantly improves the fit compared to the previous 
thermal model, resulting in \ 2 / v = 135/145, only slightly 
worse than the hybrid fit. The fitted plasma parameters are 
kT = 61 keV, r = 2.1, and kT = 83 keV, r = 1.1 for the 
two components. Indeed, plasmas in an accretion flow are 
expected, on average, to increase its temperature with de- 
creasing optical depth as a result of less effective Compton 
cooling, in agreement with our fit results. Thus, the data 
set 4 is as consistent with emission from regions of differ- 
ent kT, t as with the presence of a non-thermal tail in a 
uniform electron population. The latter interpretation may 
be favoured by the similarity of our result to those of Mc- 
Connell et al. (2000, 2002) for Cyg X-l. 

On the other hand, either thermal or hybrid model 
applied to the data set 3 yields unacceptable fits, e.g., 
X 2 jv — 202/131 in the thermal case. However, both 
the OSSE lightcurve and that from the All-Sky Moni- 
tor on board RXTE (available on the public database at 
xte.mit.edu/ASMJc.html) have shown that the source was 
significantly variable on timescale of days during the ob- 
servation. In particular, the OSSE spectrum underwent a 
significant softening and brightening during the two-week 
observation (Smith et al. 1999). The softening implies ei- 
ther kT or r of the plasma decreasing with time, to that 
during the RXTE observation (three days after the OSSE 
one). We have accounted for this effect by allowing for dif- 
ferent values of r and the same kT (since the RXTE data 
alone constrain it only weakly) in the fits to the RXTE and 
the OSSE data. We find the best-fit r for the OSSE data of 
2.4, indeed larger than the value for the RXTE data (Ta- 
ble |^). Then, the model becomes statistically acceptable, 
with x 2 / u ~ 159/130. The relative normalization of the 
OSSE and PCA data corresponds to the ratio of the respec- 
tive model luminosities of 0.67, while the normalizations of 
the HEXTE A and B data, relative to the PCA data, were 
0.891^4 and 0.82±°;°I[, respectively. 

Our results differ significantly from those of Smith et al. 
(1999) and Bottcher et al. (1999). They analyzed the same 
data and concluded that the RXTE-OSSE spectrum could 
be modelled without reflection. In order to reproduce the 
hardening at ~ 10 keV (characteristic to reflection), they 
introduced two different power-law models intersecting at 
~ 10 keV. We consider this approach (analogous to that 
of Dove et al. 1998 applied to Cyg X-l) to represent an 
only phenomenological (i.e., unphysical) description of the 
spectrum. The hardening at ~ 10 keV is, at very high prob- 
ability, indeed due to reflection (and not due to a fortitous 
intersection of two power law components just at 10 keV), as 
it is present in many other spectra of GX 339-4 as well as in 
many other Galactic and extragalactic accreting black-hole 
sources (e.g., Gilfanov et al. 2000; Nandra & Pounds 1994). 
In fact, the models of Smith et al. (1999) and Bottcher et 
al. (1999) do include an Fe Ka line and an Fe K edge, which 



features naturally appear in the presence of Compton reflec- 
tion. 

In the X-ray data (from both Ginga and RXTE), the 
reflector was found to be moderately ionized, with £ ~ 10 2 
erg s~ cm (see discussion of physical implications in Z98). 
In the Ginga data, E-p e was > 6.4 keV and consistent with 
modest ionization. On the other hand, it was < 6.4 keV in 
the RXTE data, where disc ionization was also high. Such 
an effect was also found in other RXTE observations of GX 
339-4 (Wilms et al. 1999; Feng et al. 2001), which may indi- 
cate a calibration problem. The equivalent width of the line, 
Wf c ~ 70 eV or so, is roughly consistent with the observed 
reflection (George & Fabian 1991; Zycki & Czerny 1994). 
We also see in Table [] that relativistic smearing of the re- 
processing features cannot be constrained by the data. Only 
the RXTE data appear to constrain r- ln to < 10 2 . 

4.3 A possible relation between the high-energy 
cutoff and luminosity 

The four bright-state observations of GX 339-4 cover a fac- 
tor of ~ 3 in the luminosity, L. On the other hand, all the 
spectra have very similar values of the X-ray slope, F ~ 1.75, 
as seen in Table 0. An interesting issue is then a possible cor- 
relation between L and the high-energy cutoff in the spec- 
trum. Such a correlation is generally predicted by theories of 
accretion flows, as a higher accretion rate, M, usually yields 
higher both L and r. If the X-ray slope is constant, the 
electron temperature (determining the cutoff) should then 
decrease, an effect due to an increased efficiency of cooling 
per electron. 

A complication for our data arises due to the possible 
presence of non-thermal electrons. A purely thermal fit to 
a hybrid plasma will overestimate the actual kT, as can be 
seen by comparing the results of the fits 4a and 4b in Table ^ 
On the other hand, both thermal and non-thermal electrons 
Compton scatter photons, and the average electron energy, 
(E), is then the crucial factor determining the low-energy 
slope of the Comptonization spectrum and the position of 
the high-energy cutoff. These effects are taken into account 
in the COMPPS model. We can then look directly for a cor- 
relation between L and r using fits with that model. Fig. 

shows that r indeed increases almost monotonically with 
increasing L provided the hybrid fit is used for the data 
set 4, supporting the presence of a non-thermal tail in that 
observation. Similarly, (E) shows a decreasing trend with 
increasing L, as shown in Fig. ^b. (See also Zdziarski 2000 
for a comparison of the spectra from the sets 1 and 4.) 

Zdziarski (1998) considered a simple hot accretion 
model parametrized by the Compton parameter, y (= 
ATkT/m e c? in the thermal case, and ~ 1 for our data). The 
Compton parameter determines in turn the X-ray slope, Y, 
which is constant among the 4 data sets. For constant y, 
r oc L 2//7 and L 1 ^ 6 in the advection and cooling dominated 
cases, respectively (Zdziarski 1998). The two dependences 
are shown in Fig. ^a by the solid and dashed curves, respec- 
tively. We see that the former appears to fit the data slightly 
better, but the result is not conclusive. The corresponding 
relations for (E) are shown in Fig. ^b, where we assumed 
V oc (E)t. 

Physically, a constant y (or T) corresponds to a con- 
stant ratio of the Comptonization luminosity to that in the 
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Figure 3. Relation between the luminosity and (a) the Thomson 
optical depth, and (b) the average electron energy. Labels corre- 
spond to the fit numbers from Table |^, and the dotted line denotes 
the purely thermal model to the data set 4, which is statistically 
less likely than the hybrid model. Solid and dashed curves show 
the fitted dependences on L in the cases of ady ection and cooling 
dominated flows, respectively (see Section 4.3). 



seed photons (e.g., Beloborodov 1999b), which can be due 
to a constant geometry of the flow. We also note that the 
values of L in our data are close to the maximum in a hot 
accretion flow with advection and Comptonization cooling, 
L w 10 3 V /5 (a v /0.1) 7/5 (Af/10M o ) (Zdziarski 1998), where 
q v is the viscosity parameter. 



5 THE OFF-STATE SPECTRUM 

During the fifth observation (Table |l|), GX 339-4 was by 
2 orders of magnitude fainter than during the hard-state 
observations. At this flux level, GX 339-4 has been classified 
as being in the off state, see, e.g., Kong et al. (2000). We 
note here that the object also enters even more quiescent 
states, e.g., during an observation of Asai et al. (1998), when 
an upper limit on the source flux was about two orders of 
magnitude below the flux of our data. 

We first fit the data with the same thermal model as 
for the sets 3-4, in which we assume the relative normaliza- 
tion of the OSSE data to the PCA ones of unity. This pro- 
vides a statistically acceptable model, with x 2 l v — 46/64, 
unconstrained reflection, and a broad range of the allowed 
electron temperature, kT — 1701 23 q keV. However, the spec- 
trum shows a distinct hardening at ~ 7 keV, which in turn 
requires either a very high temperature of the seed pho- 
tons, kT scc d ~ 1 keV, or very strong Compton reflection, 
S1/27T 3> 1. The former is rather unlikely at the observed 
low luminosity level, L — 8 x 10 34 erg s _1 , and the latter 
has never been observed in any other black-hole binary. We 
have performed extensive tests checking also the possibility 
that the hardening is due to the PCA background subtrac- 
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Figure 4. The upper panel shows the off-state spectrum (er- 
ror bars) from the fifth observation, modelled by thermal 
bremsstrahlung (dashed curve), power-law Comptonization at low 
energies (dotted curve) , and an Fe K line (dot-dashed curve) . The 
solid curve gives the sum. Lower panels give the contribution from 
data bins to the total \ 2 ana the data to model ratio. 



tion uncertainties. We have concluded it to be unlikely and 
hereafter consider the hardening to be intrinsic to the source, 
but still cannot rule out completely a significant effect of the 
background subtraction. 

On the other hand, it is possible that spectrum be- 
low several keV is the standard Comptonization power law, 
and the hardening at higher energies is due to thermal 
bremsstrahlung. That process is indeed expected to domi- 
nate the high-energy spectrum at very sub-Eddington accre- 
tion rates, see, e.g., Narayan & Yi (1995), Esin, McClintock 
& Narayan (1997), Zdziarski (1998). We model the spectrum 
by the sum of a power law (corresponding to Comptoniza- 
tion at low energies), thermal bremsstrahlung and a Gaus- 
sian Fe Kq line. This yields a good fit, with x 2 l v = 46/66. 
The slope of the low-energy power law is V — 3.1^'^ whereas 
the bremsstrahlung temperature is unconstrained. We have 
assumed it to equal 200 keV (e.g., Esin et al. 1997), which 
corresponds to L ~ 10 35 erg s _1 , and, e.g., r ~ 0.35 and 
the source radius of R w 10 9 , with the bremsstrahlung lu- 
minosity scaling approximately as T x ^ 2 t 2 R. The data and 
the model spectrum are shown in Fig. [jj. 

An important issue here is the origin of the strong Fe 
Ka line seen in the data. In the model with bremsstrahlung, 



the line is narrow at E Fe = 6.8+q' 3 keV with the photon flux 
of 3.2 x 10" 5 cm" 2 s" 1 , corresponding to W Fe = 530+260 
eV. If the line were intrinsic, its energy indicating origin in 
a strongly ionized medium and the large equivalent width 
would be difficult to explain given the very low luminosity 
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of the source, when the accretion disc is expected to be cold 
and far away from the black hole. For example, the Fe K 
emission is weak in a faint hard state of the black-hole binary 
XTE J1118+480, at an Edddington ratio somewhat higher 
than that of the present observation (Frontera et al. 2001b). 

Thus, we have investigated whether the line may be not 
intrinsic to GX 339-4. In fact, the source lies close to the 
direction of the Galactic centre, where the Galactic diffuse 
emission is strong. Based on the observations of Yamauchi 
& Koyama (1993), we have estimated the intensity of the Fe 
6.7 keV line in the Galactic diffuse emission at the direction 
of GX 339-4 is» 3xl0 -5 cm" 2 s" 1 deg -2 . Given the field of 
view of the PCA of ~ 1 deg 2 , the Fe K line seen in our data 
should be attributed mostly to the Galactic background. 

Similarly strong lines were reported and attributed to 
GX 339-4 in a number of off-state PCA observations (includ- 
ing ours) by Feng et al. (2001), but we believe their origin is 
also diffuse. Furthermore, Nowak et al. (2002) found that the 
ratio of the line equivalent width to the relative reflection 
strength in GX 339-4 was increasing with the decreasing 
luminosity, which they claimed to represent an argument 
against simple reprocessing models of those features. 

However, we have checked that the effect seen by Nowak 
et al. (2002) can be explained by the relative contribution of 
the diffuse 6.7 keV line increasing with the decreasing lumi- 
nosity. For the two weakest of their observations, where the 
3 — 9 keV source flux was by factors of ~ 13 and ~ 40 higher 
than in our observation, the total line equivalent width was 
Wfc ~ 150 eV and Wfb = 160 eV, respectively. The Galactic 
diffuse emission would then contribute to the observed Wfc 
at the level of ~ 40 eV and w 13 eV, respectively. Then, 
the agreement between the intrinsic Fe Ka line equivalent 
width and the strength of reflection in their data would be 
significantly improved. 



6 THE ORIGIN OF SEED PHOTONS AND 
CONSTRAINTS ON MAGNETIC FIELD 

Here, we explore some consequences of the possible pres- 
ence of a non-thermal electron tail, found statistically nec- 
essary in our uniform-plasma model for the fourth observa- 
tion. A weak non-thermal tail beyond a Maxwellian distri- 
bution leads to a moderate hardening of the X-ray power 
law in the Comptonization spectrum (as repeated scatter- 
ing off thermal electrons still dominates at those energies) 
as well as it yields a weak high-energy photon tail at ener- 
gies 3> kT (from scattering off non-thermal electrons), see, 
e.g., WZ01. On the other hand, it can amplify the cyclo- 
synchrotron emission by a very large factor, due to the fact 
that the unabsorbed synchrotron radiation is emitted just 
by the high-energy tail of the electron distribution, which 
itself is strongly amplified by the tail with respect to the 
Maxwellian (WZ01). This amplified synchrotron emission 
within the Comptonizing plasma provides then a copious 
source of seed photons. For given magnetic field strength, 
B, and the plasma parameters, kT and r, we can calculate 
the expected self-Compton luminosity, which can exceed the 
observed L in general. 

This provides a way to constrain B, with the maximum 
one corresponding to all the seed photons being synchrotron 
(see section 5.3 of WZ01). Since we know that a cold medium 



is also present in the plasma vicinity, which medium is also 
a copious source of soft blackbody photons, the actual value 
of B has to be even lower. 

First, we consider a hot accretion flow with nearly-virial 
ions (e.g., Abramowicz et al. 1995; Narayan & Yi 1995). Fol- 
lowing WZ00, we estimate the equipartition field strength 
within a radius of 20GM/c 2 (where most of the gravitational 
energy is dissipated) as B ~ 1.4 x 10 7 G. Then, following the 
calculations in WZ01 and using the hybrid el ectr on distribu- 
tion fitted to the data set 4 (S ~ 0.3, Section 4.2), we obtain 



the synchrotron-self-Compton luminosity of 4.3 x 10 erg 
s" 1 , which exceeds the observed one by a factor of ~ 40. The 
condition that the predicted luminosity does not exceed the 
observed L leads to B < 2.2 x 10 6 G. Although this is a 
factor of several below the equipartition, various processes 
dissipating magnetic field can lead to this condition to be 
satisfied in an actual accretion flow. On the other hand, the 
self-Compton luminosity from the purely thermal electron 
distribution fitted to the data set 4 (Section 1.2) would be 



only 1.6 x 10 35 erg s _1 , i.e., negligibly small compared to the 
observed L. Thus, the non-thermal tail amplifies thus the 
synchrotron emission by a rather large factor of ~ 3 x 10 3 . 
Consequently, most of the seed photons for Comptonization 
have to be emitted by a cold medium if the tail is very weak 
or absent. 

We have also considered the configuration of small, ac- 
tive regions above a cold disc, in which regions dissipation of 
magnetic field occurs (see WZ00 for details of the assumed 
model). Assuming that magnetic field dissipates at 10 per 
cent of the Alfven speed in 10 active regions of the size 10 6 
cm, the magnetic field required to dissipate enough energy 
to reproduce the observed L would be B ~ 10 s G. Then, 
however, the synchrotron-self-Compton luminosity from the 
hybrid plasma would be ~ 20 times the observed L. It would 
be < L for B < 2.1 x 10 7 G, which field, however, would not 
be able to dissipate enough energy to power the observed L. 
We thus conclude that the non-thermal electron tail fitted 
to the data rules out the active regions model for the hard 
state of GX 339-4. The same conclusion was obtained by 
WZ01 for the hard state of Cyg X-l. Also, Z98 found the 
model of active regions for the hard state of GX 339-4 to be 
ruled out based on independent spectral considerations. 



7 THE ORIGIN OF THE OPTICAL EMISSION 
IN THE HARD STATE 

GX 339-4 in an X-ray hard state was observed by the Ariel 
6 satellite (Ricketts 1983), which observation was in part 
simultaneous with optical observations, which found the 
source in a very bright state (Motch, Ilovaisky & Chevalier 
1981, 1982). The simultaneous data from 1981 May 28 are 
analyzed by M83. The 1-50 keV X-ray spectrum of Ricketts 
(1983) is, in fact, very similar to our Gmga-OSSE spectra. 
Given our considerations of the origin of the seed photons 
for Comptonization in the X-ray regime, it is of great in- 
terest to check whether the optical emission can be due to 
the cyclo-synchrotron process in the X-ray emitting plasma. 
Such origin of the optical emission was suggested by Fabian 
et al. (1982), but those authors have not performed any cal- 
culations of the actual efficiency of that process. 

Here, we assume that the plasma temperature and op- 
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Figure 5. The dereddened IR/optical spectrum (circles) of GX 
339-4 from 1981 May 24-28 (circles) and the 1-50 keV power-law 
fit (thick line) to the 1981 May 30-31 observation. The thin solid 
curve shows the thermal Comptonization and reflection compo- 
nents of the model to our data set 1 (from Ginga and OSSE), and 
the dashed curve shows the Comptonization component alone, in- 
cluding its extrapolation to low energies. The thin dotted curve 
shows a self-absorbed thermal synchrotron spectrum. 



tical depth during the 1981 observation were the same as 
those obtained from our X-ray data and consider in de- 
tail the requirements on other parameters implied by the 
cyclo-synchrotron origin of the optical emission. Fig. ^Jshows 
the unabsorbed spectrum of GX 339-4 from the IR to soft 
7-rays. The IR/optical data are from Motch et al. (1981), 
dereddened assuming E(B — V) = 1.2, and the X-ray spec- 
trum is from Ricketts (1983). 

Fig. |i] also shows the Ginga-OSSE spectrum. We see 
that the spectra from Ginga-OSSE and Ariel 6 are indeed 
quite similar. Furthermore, the extrapolation to the optical 
range of the thermal- Compton spectrum fitted to the Ginga- 
OSSE spectrum lies just a factor of several below the peak of 
the optical emission. Thus, the optical emission can indeed 
provide enough seed photons for Comptonization in the X- 
ray emitting plasma. 

We then assume that the uppermost point in the optical 
spectrum, EF E ~ 9 x IO" 10 erg cm" 2 s" 1 at E ~ 2.7 eV, 
corresponds to cyclo-synchrotron emission at the turnover 
frequency (i.e., where the plasma becomes optically thin to 
synchrotron self-absorption) . The corresponding spectrum is 
shown by the dotted curve in Fig. |B|. In the purely thermal 
case and with the plasma parameters as in Table ^, fit 1, the 
required magnetic field strength and the size of the plasma 
(assuming spherical geometry, and calculated as in WZOO) 
are B ~ 10 7 G and R ~ 10 9 cm, respectively. (These values 
are each an order of magnitute higher than those given by 
Fabian et al. 1982.) The required source size is rather large, 
~ 10 3 GM/c 2 . It can be reduced if either the magnetic field 
were larger (corresponding to the synchrotron turnover peak 
at both higher photon energy and flux, which is allowed by 
the data), kT were higher (as its value during the 1981 ob- 
servation is unknown), or the electron distribution contained 
a non-thermal tail. T he l ast possibility is, in fact, suggested 
by our data (Section |4.2j ), and we consider it in detail. We 
assume the electron power law index of p = 4, as before, and 
calculate the required values of B and R as a function of the 
fraction of energy in the non-thermal electrons, 8. Results 
are shown in Fig. We see that R is reduced by a factor of 
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cyclo-synchrotron emission parametrized by the relative fraction 
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several already for S ~ 5 x 10~ 4 , and for 5 ~ 0.1, both R and 
B are reduced by an order of magnitude. The field expected 
in a hot accretion disc flow in equipartition with nearly- virial 
ions would be B ~ 10 6 G and ~ 10 5 G at at R ~ 10 s and 
10 9 cm, respectively. Thus, R ~ 10 2 GM/c 2 ~ 10 8 cm, the 
equipartition B ~ 10 6 G, and 5 ~ 0.1 represent our pre- 
ferred set of the plasma parameters of the flow explaining 
the optical emission. 

Interestingly, M83 found an anticorrelation between the 
optical emission and the 1-13 keV X-rays. An anticorrelation 
can appear, e.g., if the cyclo-sychrotron emission is variable 
(due to variability of some plasma parameters) but the total 
luminosity is approximately constant. The spectral variabil- 
ity would then show a pivot between the optical and X-ray 
ranges, which naturally leads to an anticorrelation between 
the respective emission. Furthermore, M83 found that the 
optical emission leads the X-rays by 2.8 ± 1.6 s. At kT = 46 
keV (Table ^), it takes ~ 30 scatterings to upscatter a 2.7 
eV photon to 10 keV. At r ~ 3, the source size required 
for the time lag to be due to the delay during Compton 
upscattering is then ~ 10 9 cm, still a plausible value. M83 
also found a QPO at ~ 20 s in their data, which, interest- 
ingly, corresponds to the Keplerian frequency at a radius of 
2 x 10 9 cm (at M = IOMq). A somewhat lower QPO pe- 
riod, ~ 3 s, was found in the data set 4 (R01). The time 
scale of flares of 20 ms found by M83 can then correspond 
to the light travel time across the source. A remaining un- 
resolved issue would be then the lack of an anticorrelation 
found by M83 between the optical photons and those in the 
13-20 keV range. Possibly, that energy range is dominated 
by Compton reflection showing a different timing behaviour 
than the primary continuum. 

On the other hand, Fabian et al. (1982) proposed that 
the anticorrelation is due to the X-rays below 13 keV being 
due to emission of colder clouds surrounding the hot plasma 
(emitting both the optical emission and X-rays above 13 
keV). However, present fits to X-ray spectra of GX 339-4 
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as well as other black-hole binaries in the hard state tend 
to rule out this possibility, indicating that the entire X-ray 
emission is dominated by a single power law (e.g., U94; Z98, 
Gierlihski et al. 1997). 

We note that similar optical/X-ray behaviour has been 
observed from the black-hole binary XTE J1118+480 (Kan- 
bach et al. 2001). Optical emission from that object is very 
strong, comparable to X-rays, similar to the case of GX 339- 
4. Also, a similar anti-correlation between optical and X-ray 
flux was observed, with the optical minimum leading the 
X-rays maximum by ~ 2 s. Furthermore, the optical emis- 
sion has a narrower auto-correlation function than that of 
X-rays, which implies that the primary variability is in the 
optical band. It is then possible that the optical emission 
is cyclo-synchrotron, which is then Comptonized up to the 
X-ray range. The anticorrelation can be explained as in GX 
339-4, by spectral pivoting between the optical and X-ray 
bands. 

The site of the hot plasma could be in principle an out- 
flow, possibly a jet, which presence is indicated by radio ob- 
servations in the hard state (e.g., Fender et al. 1997, 1999; 
Corbel et al. 2000). In fact, double-peaked emission lines 
have been found in the soft state of GX 339-4 by Soria et 
al. (1999), whereas they observed broad, single-peaked pro- 
files in the hard state. They attributed the difference in the 
line profiles to the origin of the line in the outflow in the hard 
state and in the disc in the soft state. Also Wu et al. (2001) 
reached a similar conclusion. On the other hand, there is no 
natural interpretation to the frequency of the optical QPO 
in the outflow model. A way to test the origin of the optical 
radiation would be to measure its polarization, which would 
be non-zero in an outflow with a directed magnetic field. 

While the optical spectrum may in principle be either 
quasi-thermal and produced in an accretion flow or purely 
non-thermal and produced in a jet, our fits to the 7-ray 
high-energy cutoff confirm that the distribution of electrons 
producing the main X/7-ray continuum in the hard state 
has a predominantly thermal character (even if it does orig- 
inate in an outflow, as, e.g., in the coronal outflow model 
of Beloborodov 1999a). This represents a strong argument 
against the model of the X/7-ray continuum of black- hole 
binaries as being predominantly non-thermal (Markoff, Fal- 
cke & Fender 2001). In that model, the predicted position 
of the cutoff in the photon spectrum is very sensitive to 
the magnetic field strength and details of acceleration and 
cooling. Markoff et al. (2001) deal with this problem by fine- 
tuning the high-energy cutoff in the distribution of the non- 
thermal electrons. However, virtually all black-hole binaries 
in the hard state have remarkably uniform high-energy cut- 
offs (e.g., Grove et al. 1998; Z98; Gierlihski et al. 1997), 
which thus represents a severe problem for the non-thermal 
model. On the other hand, this uniformity is accounted for in 
the thermal (or quasi-thermal) model, in which the Comp- 
ton cooling (e.g., Zdziarski 1998) and/or e pair production 
(e.g., Malzac, Beloborodov & Poutanen 2001) naturally con- 
strain the electron temperature to ~ 50-100 keV. 



8 X-RAY VARIABILITY 

The fractional variability of GX 339-4 in X-rays was in- 
vestigated by, e.g., Maejima et al. (1984), Nowak et al. 
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Figure 7. The fractional rms X-ray variability of GX 339—4 in 
the PCA data of the set 4 for four ranges of frequency. 



(1999), Lin et al. (2000) and R01 and was found to be, in 
general, energy-dependent. Here, we investigate the energy- 
dependence of the X-ray variability in the RXTE data of the 
set 4, which has the best statistics of the data sets studied by 
us. We use all layers in all available PCA detectors (PCUs 
0-4) in order to obtain the best possible photon statistics. 

We construct light curves in a number of energy bands 
with different bin size, AT. For each lightcurve, we inte- 
grate the power density spectrum in the frequency range 
from 0.001 Hz up to the upper frequency of 1/(2 AT), from 
which we obtain the (Poisson-corrected) rms variability di- 
vided by the mean count rate. The latter has been corrected 
for the instrumental background, which is strong,especially 
at E > 20 keV. The results are presented in Fig. fa 

A common feature of all the shown dependencies is a 
decrease of the rms/mean starting at E ~ 7 keV. This agrees 
with a general trend of the fractional rms decreasing with 
increasing X-ray energy in other black-hole binaries in the 
hard state (Revnivtsev et al. 2000) . In our case, the ratio of 
the rms at E S> 7 keV to that at E <C 7 keV increases with 
the increasing upper frequency. 

R01 have studied a number of PCA data for GX 339-4 
including our set 4 and found that the reflection component 
is much more variable at <0.1 Hz than at ~ 10 Hz. This 
effect will result in a reduction of the rms at _E>10 keV. 
However, we have checked that the dependencies of rms on 
E in Fig. Q cannot be explained just by a superposition of 
different values of the rms for the power law and for re- 
flection. Also, the ratio of the rms/mean at high energies 
to that at low energies would decrease with the increasing 
upper frequency whereas we observe the opposite behaviour. 

It is thus more likely that we observe a change of the 
variability pattern in the Comptonization continuum itself. 
We propose that it could take place in a similar scenario 
to that considered in Section ^ when analyzing the optical 
and X-ray variability. Namely, if the total Comptonization 
luminosity remains constant but the flux in seed photons 
varied, this would result in a pivoting of the spectrum. The 
pivot energy would in general depend on the spectral index 
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and the temperatures of the electrons and the seed photons. 
Because a part of the spectrum closer to the pivot energy 
would be less variable than that further away, this scenario 
can qualitatively reproduce the observed rms/mean for the 
pivot at E > 20 keV. The dependence on the upper frequency 
can be explained if the pivoting occured mostly on long time 
scales. We note that a similar pivoting is indeed observed in 
Cyg X-l on very long time scales (Zdziarski et al. 2002). 

We note that the pivoting energy required by the 
above interpretation is different from that needed for ex- 
plaining the optical/X-ray anticorrelation by varying cyclo- 
synchrotron seed photons in Section 0. For that, the pivot 
energy has to be obviously between the optical and X-ray 
ranges. Thus, the two interpretations are mutually exclusive. 
We have also calculated that the pivot at E > 20 keV im- 
plies fcTseed 0.1 keV, which rules out the cyclo-synchrotron 
origin of the seed photons in our data set 4. We point out, 
however, that the two variability patterns were observed in 
two different observations of GX 339-4 (separated by 15 
years), during which the state of the source could be differ- 
ent. Also, the cause of the dependence of the rms on E shown 
in Fig. (7] could be different than due to pivoting. It can pos- 
sibly be explained, e.g., by models of magnetic avalanches 
(e.g., Poutanen & Fabian 1999), which predict complicated 
patterns of variability of flares on different time scales. 



9 CONCLUSIONS 

We have analyzed four broad band X/7-ray spectra of GX 
339-4 in the hard state and found that they are well fitted 
by Comptonization being predominantly thermal, Compton 
reflection and fluorescent Fe K emission with the strength 
compatible with that of the reflection. This provides a fur- 
ther confirmation of the dominance of those physical pro- 
cesses in this object, and black-hole binaries in the hard state 
in general. It also argues against the non-thermal model for 
the origin of the X/7-ray continuum of black-hole binaries 
by Markoff et al. (2001). We also study time variability in 
our best X-ray data set, and find a decline of the fractional 
rms with increasing photon energy, again consistent with 
previous findings for this source. 

Our fifth observation was in an off state, with the X- 
ray flux two orders of magntitude below that in the hard 
state. Our major finding here is that the very strong Fe 
K line present in our data (as well as in other reported 
proportional-counter off-state observations of this source) is 
not intrinsic to GX 339-4, but is due to the strong diffuse 
6.7 keV emission in the Galactic Center region. 

We also obtain a number of other results, which, how- 
ever, are less unambiguous and require future further stud- 
ies. This is due to either the limitations of the available data 
or insufficient theoretical understanding at the present time. 

In the off state, we find a spectral hardening at ~ 7 
keV, which we find to be incompatible with thermal Comp- 
tonization but, rather, compatible with the presence of a 
bremsstrahlung spectral component. Such a component is in 
agreement with theoretical predictions, but its unambiguous 
detection is hampered by the low signal-to-noise ratio and 
lack of sufficient coverage at high energies. 

In the hard state, we find a statistically significant de- 
parture from the pure thermal Comptonization spectrum at 



high energies in our best spectrum. We interpret it as due 
to the presence of a weak high-energy tail in the electron 
distribution. However, we cannot rule out an alternative ex- 
planation of the effect of the plasma parameters varying in 
space and/or time. 

The possible presence of a high-energy electron tail has 
implications for the strength of the magnetic field in the 
source, from the consideration of Compton scattering of the 
cyclo-synchrotron emission. The magnetic field strength is 
constrained to values somewhat below equipartition (with 
hot ions) in the case of a hot accretion disc. 

Interestingly, the four hard-state spectra have very sim- 
ilar intrinsic X-ray slope, V ~ 1.75. On the other hand, the 
high-energy cutoff energy decreases with increasing luminos- 
ity, and, correspondingly, the fitted Thomson optical depth 
of the plasma increases with luminosity. The significance 
of the correlations is strongly increased if we allow for the 
presence of a high-energy electron tail in our model. Such a 
correlation is in agreement with theoretical predictions for 
hot accretion flows. Still, it is based only on the four avail- 
able spectra, and studies of more broad-band spectra with 
similar X-ray slope, F, are highly desirable. 

We also study the origin of the strong optical emission 
found in the past to be associated with the X-ray bright hard 
state. We find the reported optical emission can be due to 
the cyclo-sychrotron process, and the emission also provides 
seed photons for Comptonization. The presence of a high- 
energy electron tail alleviates constraints on the emission 
region, allowing it to be smaller and with a weaker field. 
The reported optical/X-ray anticorrelation can be due to 
the synchrotron emission being variable at an approximately 
constant total luminosity, which leads to a spectral pivot 
between the optical and X-ray bands. 

On the other hand, we find the dependence of the rms 
variability on energy can be also explained by variable seed 
photons, but the pivot point is then at > 20 keV, i.e., in con- 
flict with the previous model. Thus, a comprehensive model 
of both of these variability patterns requires more compli- 
cations, probably more source components. 
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